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Nuclear effects in polarized inelastic electron scattering off polarized ^He and po- 
larized are discussed; in the resonance region, Fermi motion effects are found 
to be much larger than in deep inelastic scattering. It is shown that improperly 
describing nuclear dynamics would lead to the extraction of unreliable neutron spin 
structure functions; on the other hand side, simple and workable equations, relat- 
ing the Gerasimov — Drell — Hearn Integral for the neutron to the corresponding 
quantity for ^He and ^H, are proposed. Nuclear effects in the recent E143 data 
are estimated by a proper procedure. 

1 Introduction 

The measurement of the polarized nucleoli Spin Structure Functions (SSF) gi 
and §2 in the resonance region allows one to check the helicity structure of the 
photon ~ nucleon coupling between the Deep Inelastic Scattering (DIS) region 
and the real photon limit El. Recently it has been proposed at TJNAF to study 
the SSF in a wide range of energy (0.2 GeV < v < i Ge\Q and momentum 
(0.15 GeV'^ < < 2 GeV'^) transfers, for both the protoia and the neutron, 
using in the latter case polarized and '^He targetflu. 

For a hadronic target A, the SSF gi{v,Q^) and 92 {^tQ^) are exneri- 
mentally obtained by measuring longitudinal and transverse asymmetriefl. A 
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relevant quantity related to the SSF's 51(2) is the following integral 



K 



(1) 



where vth = (Q^ + 2m-^m + rn^)/2m is the threshold energy for the pion-ele- 
ctroproduction off the the nucleon, m is the nucleon mass, m-^ is the pion mass, 
K is the photon flux, a is the fine structure constant and Gi{v^ Q"^) reads as: 



I 



(2) 



Two features of Qii^^ Q^) have to be mentioned: (i) when multiplied by the 
quantity —jf^i^ + ^) , Ma being the target mass, it coincides with the usual 
transverse response, Rt' (see e.gi); and (ii) in the DIS limit {Q^ 00, v 00, 
Q^/iy fixed), coincides with the SSF, g^{iy, Q^). 

For any spin ^ hadronic target, the SSF and read as foUowfl 
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(4) 



where ^o'^iy, Q'^) — a^i^iv, Q'^) — (Ty^{v, Q^), ('■fy2(3/2) cross sec- 

tion for photon - hadron scattering with total helicity 1/2 (3/2), cr^(z/, Q^) is 
the transverse - longitudinal interference cross section. Thus Eq. (||) becomes: 

MaK 



8^2^(1 + 1^) 



^f/2(^: Q)~(^K/2{^^ Q) 



(5) 



For a nucleon target (A = iV, N = n 01 p), the integral coincides with 
the Gerasimov - Drell - Hearn (GDH)Jntegral, loDHiQ^)^ which, in the real 
photon limit, gives the GDH Sum Rulfl: 

00 
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iSdh{Q' = 0) = / ^(af/2(-,O^ = 0)-a3^/2(-,Q^ = 0)) 



-0.53 GeV ^ for protons 
-0.60 GeV'^^ for neutrons 

where k^t is the anomalous magnetic moment of the nucleon. 

An important observation^ can be made about the evolution of Iqdu- 
Since in the large limit Qi coincides with , one can evaluate (||) for the 
proton, using the resultsS at ~ 10 GeV"^, which gives IcDniQ'^) - 0.14/Q^ 
This result, if compared with (|^), provides evidence of a sharp modification 
in the helicity structure of the -fp coupling between the real photon limit 
and the DIS region, which lead to a change of sign of Iquh^Q^) some value 
of . Iii_tpdcr to understand this behavior, much theoretical work has been 
produceiyO, both in the real photon limit and at finite values of Q^. The 
electroexcitation of nucleon resonances is evidently the main reason for this 
evolution of the integral, Igdh{Q^)- Therefore, experimental investigation 
of the low-Q'^ evolution of Iqjj^ is of great relevance, in particular for the 
neutron, for which several analyses of the available unpolarized data for pho- 
toproduction and low-Q^ electroproduction disagree with the prediction (^, 
while similar estimates seem to give the correct value for the proton. 

Since Gi{v, Q^) will be extracted from experimental data on "^H and ^iJe, 
corrections due to nuclear effects have to be introduced. The goal of this talk 
is to illustrate the relevance of these corrections in the resonance region. Our 
results for '^He and "^H are presented in the following two sectionsSEj. The 
recent deuteron data of the_E143 coUaboratioiEJ are used as an example to 
test the method proposed iiEa for the extraction of the neutron SSF from the 
nuclear ones. 



2 The process ■^He{e, e')X in the resonance region 

Using the convolution model for the nuclear SSF 9i(2) described ii^'Ei0, we 
obtaiiS 

ofi-^Q') - J2 [ [ [ dp^^ {gn-'^Q') P\^{P.E) 

AT \J U U jD 1 i 

N—p,n ^ 

+ T''ip,E,Q^)}s(.+ ^^^^p^-l^) , (7) 
J \ Zmv mv J 

where E is the removal energy of a nucleon with momentum p, p = {Ma — 
y'' {E + Ma — rriY + p) is the nucleon 4-momentum, q is the 4-momentum 

transfer, Ep = y/m"^ + \p\'^, {p, E) and T (p, E, Q^) are related to 
the elerojsnts of the 2x2 matrix, which represents the spin dependent spectral 
functioiQ. 
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Figure 1: qI""" in DIS (Q^ = lOGeV^) (a)Q, and in resonance (Q^ = i GeV^) (b) regions, 
obtained by considering Fermi motion and binding (full). The dashed curve represents the 
same functions obtained considering the proton and neutron effcijtive polarization in "^He as 
the only relevant nuclear effects (Eq. (|l2[)). 



The elements of this matrix are 



/a-1 

X 5{E-E{_^+Ea) (8) 

where \iI^jm) is the ground state of the polarized target nucleus, is the 

eigenstate of the (A-1) nucleon system with eigenvalue Ej^_^, and \p, a) is the 
plane wave describing the struck nucleon in continuum. The relevant quantity 

7-1 

is r« {p, E)^ Pfi j_ (p, E) ~ P^i _ 1 1 (p, E), i. e. the difference between the 

' 2 2 '2 2 2 ' 2 

spectral functions with nucleon spin parallel or antiparallel to the nucleus spin. 
It provides the effective nucleon polarizations pjv produced by the S' and D 
waves in the ground state of ^He: 



Pn = 



dE 



N 



{P.E) 



(9) 



The spectral function om gives Pp(„) =p^0.030 (0.88), in agreement with world 
calculations on the three body systemska, pp(„) -0.028 ±0.004 (0.86 ±0.02). 

The term T (p,E,Q^) depends also upon Pi_i i{p,E), as well as upon a 

I — I 2 2 ' 2 

proper combination of the SSF's gi and g^. However, it is_of the order \p/m\, 
and thus gives a very small contribution both in the DIx3 as well as in the 
present calculation of the resonance region. For this reason it will be omitted 
hereafter. The quantity which appears in Gi^*^ (Eq. (|^)) is defined in 
terms of the polarized electroproduction cross-sections 0^2(3/2) (*^^- ^1- ®)- 
Since thejfirst few data for appeared very recently tJ, we use a theoreti- 
cal modelllJ, where the contributions of the resonances P33(1232), £)i3(1520). 
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511(1535) and i^i5(1680) have been parametrized using the existing experi- 
mental data for unpolarized electroproduction. Other resonant states have 
been added using the predictions of a single quark transition model, and the 
single pion Born term background has also been included. Using the above 

models for P|l {p,E) and , we have calculated Qi^'^ in the resonance re- 
gion. The results are presented in Fig. 1, where they are compared with 
the results for the DIS limit0. In the latter case giiiy,Q^) = gti.v.Q'^) = 

a^l^{v, Q^) — 17^/2 ('^' Q^) ' l^'-^- (0) reduces to the well-known con- 



MaK 

volution formula: 



Ma jm 



gn^,Q') = I rf^sf G^(^) , (10) 

where (z) is the light-cone momentum distribution 

G^(z) = JdE JdpP^^ip,E) s(^z-^^ (11) 

with p"*" = (p^v — p ■ q)/\q\ being the light-cone momentum component. Fig. 1 
shows that in the DIS case the following equation 

gl"^x,Q') « 2ppgpix,Q') +pngn^,Q') , (12) 

approximates very well the convolution formula, at least for x < 0.80. The 
same does not hold in the resonance region, where nuclear effects turn out to 
be important, with the Fermi motion significantly broadening and damping the 
peaks associated with the excitation of the prominent resonant states. There- 



fore, Eq. (12) can be considered as a workable formula for extracting the DIS 
5i (x, Q^) from experimental data on g\^'^(x^ Q^), but in the resonance region 
it appears to be of little help. We stress that the importance of nuclear effects 
in the resonance region is a well-known feature of the unpolarized scattering 
as well. Nevertheless, we have found that the proton contribution to the nu- 
clear Gi^'^{x,Q^) is not large, and therefore ^He is a good effective polarized 
neutron target also in the resonance region. 

Let us now discuss the role of nuclear effects on the integral Eq. (|l|) . In 
Fig. 2 the integral, calculated placing Eq. (^ in Eq. (|l]), is compared with 
the following expression 

r^-(Q2) = 2ppP{Q')+pnr{Q') , (13) 
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which represents the integral within the assumption that Fermi motion 
and binding can be disregarded and that nuclear effects in polarized ■^iJe are 
only due to the effective nucleon polarizations. It can be seen that Eq. (13) 
approximates Eq. (|l|) with an accuracy better than 5%. From the same figure it 
can also be seen that, because of the efiective nucleon polarizations, the integral 
for the free neutron predicted byliil, noticibly differs from / ^^{Q"^)- Thus 
the nuclear structure effects are very large in the nuclear SSF Qiiv^Q'^) (cf. 
Fig. 1 (b)), but in the integral / ^^(Q^) these effects can be approximately 
accounted foe, by the effective nucleon polarizations (cf. Fig. 2). This result is 
understood t2l by performing a series expansion of the r.h.s. of Eq. (j^) in the 
variable z around the non relativistic value z — 1. 

To sum up, we have shown that in the resonance region Qi^''{x,Q'^) ^ 
PnGi{x,Q'^) +2ppQf{x,Q'^), but the integrals of the two quantities are very 
similar. 



3 The process ^i?(e, e')X in the resonance region 

The nucleon contribution to the deuteron structure functions is usually calcu- 
lated by weighting the amplitude of electron scattering on the nucleon with th e 
wave function of nucleon in the deuteron (for recent developments see e.g.EjEij 
and references therein). For the SSF the most important effects are the Fermi 
motion and the depolarizing effect of the D-wave. Additional effects, jSiich as 
off-mass-shell effects or nucleon deformation, are found to be small E3. For 
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finite values of and v, Eq. (|^) for the deuteron reads as followH 

Q') - I i^-G^ {^\Q') { 1 + ^^^^ 

J (2iTy pq ^ ' \ m 

-g^{x/z,Q^)fD{z,^{x,Q')), (14) 

where = {Q^ + Q^)/2 is the isoscalar nucleon response given by Eq. (||) 
and '^'^{p) the deuteron wave function with spin projection M. In the rest- 
frame of the deuteron, with q opposite the z-axis, the kinematical variables 
are defined as follows: 

pq iy{po + £,{x,Q'^)p3), po = m + ed - p'^ /2m, (15) 
^ = qs/iy = \q\/u = + 4m2a;Vg2, = V, x* = Q"" /2pq, (16) 

where eo = —2.2246 MeV is the deuteron binding energy. 

The limits Zrain{max){x ^ Q^) ^re defined to provide an integration over the 
physical region of momentum in ( p^ and to take into account the pion pro 



duction threshold in the virtual photon- virtual nucleon scattering. Eq. (14) 
has the correct limit in DIS. In this case: ^{x, Q^) 1, Zmin — > x, z^ax — > 
Mn/m, aiuLllie usual convolution formula for the deuteron SSF gi{x,Q'^) is 
recovered c3Eil: 

Mo /m 

5f(x,Q2)^ 7 ^g^{x/z,Q')fn{z). (17) 



(cf. Eq. (12) for the ^He case). Equation (17) defines the spin-dependent 
"effective distribution of the nucleons" , fo, which describes the bulk of nuclear 
effects in . The main features of the distribution function, /^{z), are a 
sharp maximum at z = 1 -I- e£)/2m « 0.999 and a normalization given by 
(1 — 3/2Pd) (Pd being the weight of the D-wave in the deuteron). As a result 
in the region of medium values of x ~ 0.2 — 0.6 the deuteron SF gP{x) is 
slightly suppressed by Fermi motion and binding effects, compared to (1 — 
3/2Pd) X g^{x). However, the magnitude of this suppression is small (~ 1%) 
and this is why it is phenonienologically acceptable to extract the neutron SF 
from the deuteron and proton data by making use of the following approximate 
formula: 

9?{x, Q') ^Pd) {gn^, Q') + gUx, Q'))/2 , (18) 
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(cf. Eq. (14) for the ^He case). In addition, when integrated over x, eqs. (17) 
and ( p^ ) give exactly the same result (F = J dxgi{x)), i.e. 

TDiQ') =(l- Ipd^ (r„(Q2) + Fp(Q2))/2, (19) 

which allows one to obtain exactly the integral of the neutron SF F„ knowing 
the dcuteron and proton integrals, without solving (pT|). 

As in the '^He case, the equations at finite values of and v arc more 
sophisticated than the corresponding equations in the deep inelastic limit. In 
particular, Eq. ( p^ ) does not represent a "convolution formula" in the usual 
sense, since the effective distribution function fo and the integration limits 
are also functions of x. This circumstance immediately leads to the conclu- 
sion that, in principle, when integrals of the SF are considered, the effective 
distribution can not be integrated out to get a factor similar to (1 — 3/2Pd) 
in ([T^). Another interesting feature of Eq. ( p4| ) is the Q^-dependence of fn 
and (ma2.)(a;, Q^). If we again limit ourselves to the discussion of the in- 
tegrals of SF, one concludes that the Q^-dependence of such an integral is 
governed by both the QCD-evolution of the nucleon SF and the kinematical 
Q^-dependence of the effective distribution of nucleons. Thus, in principle, in 
the non-asymptotic regime, equation (p^), does not hold. 

We have performed a realistic calculation of Eq. dlj). As in the ^i?e case, 
in our numerical estimates we have evaluated i]p-(|l4|) using the elementary 
cross sections for the proton and neutron giYen inEJ. Using the Bonn poten- 
tial model for the deuteron wave functioned, we have carried out a realistic 
calculation of Qi'{x,Q^), (p^, in the region of nucleon resonances. In Fig. 
3, the obtained Q['{x,Q^) at — 1.0 GeV^, is compared with the input of 
the calculation, i.e. the isoscalar nucleon response, {x,Q^). It can be seen 
that the role of nuclear effects in the resonance region is much larger (up to 
~ 50% in the maxima of the resonances!^ than in the deep inelastic regime 
(~ 7 — 9%, depending upon the models EBlEil, with resulting 6 — 7% from the 
depolarization factor (1 — 3/2Pd) and ~ 1 — 2% from the binding effects and 
Fermi motion). As in the ^He case, such a drastic effect is a consequence of 
the presence of the narrow resonance peaks in the nucleon {x, Q^). 

Let us now discuss the role of nuclear corrections in the analysis of the 
integrals of the SSF, such as the GDH Integral. 

An important observation has been made in DIS, namely the exact formula 
( p^) and the approximate formula (^) give the same result for the integral of 
the neutron structure function, g"(a;, ^ m?) (see eq. (p^)). The applicabil- 
ity of the approximate formula in DIS is based on the conservation of the norm 
of the distribution f{z) by the convolution formula (^^)- This circumstance 
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Figure 3: Gi{x,Q'^) at Q-^=l GeV'^ for the deuteron SF (solid line), compared with the 
corresponding quantity for the isoscalar nucleon (dotted line) , used as input in the calculation 

of eq. 



can not be immediately extended to the case of the resonance region, since: (i) 
the covolution is broken in eq. ([ij) and (ii) the normahzation of the function 
f{z,x,Q^) is different from the one of f{z). The integral of the distribution 
f{z,x,Q^) represents the "effective number" of nucleons "seen" by the virtual 
photon in the process when the virtual photon is absorbed by the nucleon and 
at least one pion is produced in the final state (it is less than 1 at low and 



However, the use of the formula (|19|) in the resonance region gives results 
numerically; very close to the integration of the exact equation (^^. As ex- 
plained irt^l, this is a consequence of the smallness of the effects breaking the 
convolution in eq. (p^. 

We have found that the integrals of the SF, such as the GDH Sum Rule, 
can be estimated with accuracy better than 3% by the simple formula (19) 
which is also valid in deep inelastic region. 



4 Neutron SSF from nuclear data 



At this point, we have observed that nuclear effects in the resonance region 
are very specific and therefore approximate formulae, such as ( p^ ) for '^He or 
( |l8| ) for ^iJ, do not work even for crude extraction of the neutron response. 
Obviously, another method of extracting should be used. 

In ref. E3 a rigorous procedure of solving eq. (^|) for the unknown neu- 
tron SSF has been proposed and applied in the deep inelastic region. It has 
been shown that this method, which works for both spin-independent and 
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spin-dependent SSF, in principle allows one to extract the neutron SSF ex- 
actly, requiring only the analyticity of SSF. It can also be applied by a minor 
modification to the extraction of the SSF at finite . Such a procedure can 
be applied both to and ^He targets; here only the deuteron case will be 
discussed. 

The basic idea is to replace the integral equation ( p^ by a set of linear 
algebraic equations, KGn — Gd, where X is a square matrix (depending upon 
the deuteron model), Gd is a vector containing the available Q[' data and Gn 
is a vector of unknown solutions. Changing the integration variable in (Q), 

T = x/z, we get 

GF{x,Q^)= J dTg^{r,Q')^fD{x/T,ax,Q')), (20) 

where Tmini^X^ Q ^ — xj Zmax{x^ Q ), '^maxix-, Q ) — X^aaxi^Q ^ / '^min{x ^ Q ) and 

XmaxiQ'^) is defined by the pion production threshold in virtual photon- nucleon 
scattering. Let us assume that has been measured experimentally in the 
interval {xi,X2) and that a reasonable parametrization for it can be obtained 
in this interval. Then, dividing both intervals {xi, x^) and (Tmim T„-iax) into N 
small parts, one may write: 



N 



^f(a:„g2)«^gf(f^.,g2) / ]_f^^^^/r^Q2)dr, z = 1...7V, (21) 

where f,- = Tmm+/i(j-l/2) and h = {T„iax-T.,nin)/N . Equation (|2^) is already 
explicitly of the form Gd ^ KGn, therefore usual linear algebra methods can 
be applied to solve it. 

Note that the range of variation of r is larger than the one for x. Therefore 
the experimental knowledge of in the interval (xi, X2) provides information 
about the neutron in a wider interval (for example, in deep inelastic regime 
Tmin xjl and Tmax = 1)- Howevcr, extracting information beyond the 
interval = x\ to Tmax — X2 is almost impossible in view of the structure 
of theJcernel of eq. (|20| ) and the kinematical condition of planned experimental 
datall3. We have to redefine the kerneLof eq. (|2^) to incorporate new limits 
of integration fmin = xi and Tmax = a;2llj. The experimental errorbars of the 
nuclear data can be related to errors in the extracted structure function of the 
nucleon. 

To show how this procedure works in practice, we have obtained of the 
isoscalar nucleon solving Eq. (|lj) , using in the left hand side an analytical fit of 
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0.1 




Figure 4: The analytical fit of the deuteron E143 dataLj used as input in our example of 

extraction 



one set of the recent E143 deuteron dataiJ (corresponding to 1.01 < < 1.28 
GeV^, see Fig. 4). It should be noticed that in the experimental analysis the 
difference between gi and Qi (cf. Eq. (^)) has been neglected, and therefore the 
data refer to the SSF of the deuteron. Resonance structures corresponding 
to the three main resonances have been included in obtaining the fit. However, 
no further physical assumptions have been made, and the behavior of the fitting 
curve is basically forced by the data, whose errorbars are large. Therefore, our 
analysis is only a crude estimate of the possible nuclear effects, and presently 
it should be considered as an example of how the extraction method could 
be applied. Once we have assumed that our fit describes the nuclear data 
fairly well, Eq. ( p^ ) can be solved using the method described in this section. 
The results of the calculation are shown in Figs. 5 and 6. It can be seen 
that nuclear effects are large, and that the extraction procedure allows one to 
discover features of the isoscalar nucleon SSF which, due to nuclear effects, are 
not apparent in the initial nuclear SSF. The neutron SSF could be obtained by 
subtracting the proton data from the isoscalar SSF, and provided the errorbars 
on Q^' are reduced, this method will provide us with a reliable neutron SSF. 

5 Conclusions 

We have shown that the effects of nuclear structure in the extraction of the 
neutron SSF in the resonance region are much more important than in DIS. 
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Figure 5: The extracted ispspalar nucleoli SSF (dashed hne) compared with our fit of the 
E143 dataEj (full line), used as input in our calculation. 
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Figure 6: The extracted "data" of the isoscalar nucleon SSF (open circles) compared with 
the experimental deuteron ones (filled squares). The extracted points have been slightly 
shifted, in order to distinguish their errorbars from the the experimental ones. 
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We have explained how the correct neutron SSF can be firmly extracted from 

the combined deutcron and proton data. As for the integrals of the SSF, the 
estimates are easier, and they can be carried out by using the simple procedure 
similar to the DIS case. 
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